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Abstract. We consider two mean-field like models which belong to the 
universality class of absorbing phase transitions with a conserved field. In 
both cases we derive analytically the order parameter as function of the control 
parameter and of an external field conjugated to the order parameter. This 
allows us to calculate the universal scaling function of the mean-field behavior. 
The obtained universal function is in perfect agreement with recently obtained 
numerical data of the corresponding five and six dimensional models, showing 
that four is the upper critical dimension of this particular universality class. 
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1. Introduction 

The scaling behavior of directed percolation is recognized as the paradigmatic example 
of the critical behavior of several non-equilibrium systems which exhibits a continuous 
phase transition from an active state to an absorbing non-active state (see for 
instance Jl], |^]). The widespread occurrence of such systems in physics, biology, 
as well as catalytic chemical reactions is reflected by the well known universality 
hypothesis of Janssen and Grassberger that models which exhibit a continuous phase 
transition to a single absorbing state generally belong to the universality class of 
directed percolation ||, |J] . Introducing an additional symmetry the critical behavior 
differs from directed percolation. In particular particle conservation leads to a new 
universality class of absorbing phase transitions with a conserved field as pointed out 
in ||. In this work the authors introduced two models, the conserved lattice gas 
(CLG) as well as a conserved threshold transfer process (CTTP). The latter one is 
a conserved modification of the threshold transfer process introduced in ||. Both 
models display a continuous phase transition from an active to an inactive phase. The 
density of active sites p a is the order parameter of the phase transition controlled by 
the total density of particles p, i.e., p a > if the density exceeds the critical value p c 
and zero otherwise. As usual in second order phase transitions the order parameter 
vanishes algebraically at the transition point. The corresponding order parameter 
exponent as well as the exponent of the order parameter fluctuations of the CLG are 
determined in for various dimensions. 
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The scaling behavior of the CLG model in an external field conjugated to the order 
parameter was considered recently ||. The external field is realized by movements 
of inactive particles which may be activated in this way. Thus the field creates 
active particles without violating the particle conservation. Taking into account this 
additional scaling field the order parameter obeys the scaling ansatz 

p*(5p,h) ~ \r{8p\- l 'P,h\-°'P) (1) 

with the critical exponents j3 and a, the scaling function f, the reduced control 
parameter bp = p/p c — 1, and the external field h. Choosing SpX~ x ^ = 1 one 
gets for zero fields p a ~ ^(L0)<5p^, whereas h\~ a IP — 1 leads at the critical density 
to p a ~ r(0, V)bPI° '. Except of the critical point (bp — 0,h = 0) the scaling function 
f(x, y) is smooth and analytic but it is not universal since it may depend, like the 
value of p c , on the details of the considered systems (here e.g. the lattice structure, 
the update scheme, etc.). 

A universal scaling function R can be introduced if one allows non-universal metric 
factors Cj for the scaling arguments 5p and h (see for instance Q), i.e., 

Pa (Sp,h) ~ XR(c 1 6p\' 1/l3 ,c 2 hX-^ 13 ), (2) 

and the scaling function is normed by the conditions i?(l,0) = R(Q, 1) = 1. Then the 
function R(x,y) is universal, i.e., similar to the critical exponents R(x,y) is identical 
for all models which belong to the same universality class. But the non-universal 
metric factors differ again between the models and may depend on the lattice structure, 
the used update scheme etc. 

The non-universal metric factors can be easily determined by the scaling behavior 
of the order parameter at zero field and at the critical density, respectively. Choosing 
ci Sp \~ X /P = 1 one gets for zero fields (h = 0) 

Pa (6p,0) ~ {c x 5pf (3) 

whereas c 2 h X^ a ^ = 1 leads at the critical density (Sp = 0) to 

Pa (0,h) ~ (c 2 hf/°. (4) 

In this work we derive the universal scaling function R of the mean-field solution 
of the universality class of absorbing phase transitions with a conserved field. In 
particular we consider analytically the CLG and the CTTP with particle hopping to 
randomly chosen sites on the whole lattice. This unrestricted particle hopping breaks 
long range correlations and the scaling behavior is characterized by the mean-field 
exponents (see fl0|)- Neglecting correlations it is possible to derive analytically the 
order parameter as a function of the control parameter and of the external field. The 
obtained universal function is in perfect agreement with recently obtained numerical 
data of the five and six dimensional CLG and CTTP in an external field. 



2. The conserved lattice gas 

We consider the CLG model on a chain with L sites and periodic boundary conditions. 
At the beginning one distributes randomly N — pL particles on the system where p 
denotes the particle density. A particle is called active if at least one of its two 
neighboring sites is occupied. In the original CLG model active particles jump in the 
next update step to one of their empty nearest neighbor site, selected at random Q. 
In the steady state the system is characterized by the density of active sites p a which 



Mean-field scaling function of absorbing phase transitions with a conserved field 3 



Table 1. The configuration of a CLG lattice before (C) and after (C) a particle 
hopping. Only the target lattice site where a particle hops onto and its left and 
right neighboring sites are shown. Empty sites are marked by o, inactive sites 
are marked by *, and active sites by •. An denotes the change of the number of 
active sites due to the particle hopping and p is the corresponding probability of 
the configuration C if one neglects spatial correlations. 
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depends on p. The density of inactive sites is given by p\ = p — p a and 1 — p is the 
density of empty sites. 

We introduced in a modification of the CLG model where active particles 
are moved to a randomly chosen empty lattice site which suppresses long range 
correlations. A given lattice site is active with a probability p& and with the probability 
1 — p it may be moved to an empty lattice site. Depending on the neighborhood of 
this new lattice site the number of active sites may change. For instance if both new 
neighbors of the moved particle are empty the number of active particles is reduced 
by one, An = — 1. Without correlations the corresponding probability for this process 
is p a (1 — pY . In the case that one of the new neighbors of the moved particle is 
occupied by an inactive particle (p{) and the second neighbor is empty (1 — p), the 
number of active sites is increased by one (Art = 1). The corresponding probability is 
given by p = 2p a pi (1 — p) 2 . All other possible configurations and the corresponding 
probabilities are listed in tablc|l|. 

The probabilities that the number of active particles are changed by An are given 

by 
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The expectation value of An is 

2 

E[An] = An P^ = (l-p)Pa[-l-2p a + 4 / 9- / 9 2 ]. (6) 

An=-1 

As pointed out in JToj j the average number of active sites is constant in the stationary 
state, i.e., the expectation value of An should be zero in the steady state. Using the 
constraint E[An] — one gets 

P=l V p a = V -l-2p a + 4p-p 2 = 0. (7) 

The first equation corresponds to a system were all sites are occupied (p a = 1) and 
no dynamics can take place whereas the absorbing state is reflected by the second 
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Figure 1. The order parameter of the CLG model as a function of the particle 
density p and the applied external field h [see equation (jllj)]. The thick line 
corresponds to h = whereas the thin lines correspond to h = 0.1,0.05,0.01 
(from top to bottom). 



equation. The non-trivial third equation corresponds for p a > to the active phase 
and one gets for the order parameter in leading order 



4p - P 2 - 1 
2 



{2V3-3)S P + 0{S P 2 ) (8) 



with the critical density p c = 2 — V3 [0 . Thus we have obtained the critical exponent 
= 1 as well as the non-universal metric factor c\ = 2\/3 — 3. 

In the case that an external field is applied non-active sites may be activated 
(see Q). The probability that a site is occupied and has two empty neighbors is 
p(l — p) 2 . These particles are activated with probability h, where h denotes the 
strength of the applied field. In this process the number of active sites is increased 
(An = 1) and the probability pAn=\ is modified to be 

Pa» =1 = (1 - p) Pa [2(1 - p) Pi + 2pip a ] + (1 - pf p h. (9) 

Using again the steady state condition i?[An] = one gets the equations 

p=l V Pa [-1 -2 Pa + 4p-p 2 ] + (l-p)ph = 0. (10) 

The first equation corresponds again to the trivial case of a totally occupied lattice 
whereas the second equation yields the solutions 

1 



Pa = i (-1 + 4p - p 2 ± y/8h(l - P )p + (-1 + 4p - p 2 ) 2 ] . (11) 

The solution with the + sign describes the order parameter p a (p, h) as a function of the 
density and of the external field whereas the — sign solution yields negative densities 
for the order parameter for all values of p and h. A sketch of the order parameter for 
various fields is presented in figure |l|. 
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At the critical density p c = 2 — y/3 the order parameter is given by 



P*(Pc,h) = ^ (12) 

i.e., the field scaling exponent is a — 2 and the non-universal metric factor is 
c 2 = (3a/3-5)/2. 

In the following we derive the universal scaling function R(x, y) of the mean-field 
solution. Therefore we write the order parameter [equation (|ll|)l as a function of the 
reduced control parameter Sp and consider the function p & {5p,h) / ' \J~h. Since we are 
interested in the scaling behavior in the vicinity of the critical point we perform the 
limits p a — > 0, Sp — > and h — » with the constraint that p a /Vh and Spj\fh are 
finite. Thus all terms which scales as 5p 2 /vh or Sp\fh vanishes and in leading order 
and we get 



3^3-5 /2V3- 



p a (6p,h) 2V3 — 3 Sp 

Vh = 2 TTi + \ 

Introducing the non-universal metric factors c\ = 2^/3 — 3 and c-2 = (3-\/3 
gets the universal function 




%M = f + ^ 2 H(f) 2 . (14) 

Equations (p|,|l2|) are recovered from this result by setting h = and Sp = 0, 
respectively. Furthermore we get i?(l,0) = R(0, 1) = 1 as required above. 

As usual in scaling analysis (see for instance ||) the order parameter as well 
as the control parameter are rescaled by the field in order to obtain a data collapse 
[setting C2 h \~ a IP — 1 in equation (|J)]. In this case one gets the universal function 



Pa(<?P, h) 5, , x I fx\ 2 

where the scaling argument is given by x = c\Spj ' \[cvji. 

For the sake of simplicity we derived the scaling function of the one-dimensional 
CLG model only. A straight forward extension to higher dimensional systems for 
h = was already presented in pX| |. The increased number of nearest neighbours in 
higher dimensions affects the non-universal quantities p c , ci , and ci only, but not the 
critical exponents and the universal scaling function. 



3. The conserved threshold transfer process 

A similar analysis can be performed for the CTTP with random neighbor hopping. 
In the CTTP lattice sites may be empty, occupied or double occupied. Double 
occupied lattice sites are considered as active and one tries to transfer both particles 
of each active site to randomly chosen lattice sites. Recently performed numerical 
investigations in dimensions d — 2, 3, 4, 5, 6 confirm the conjecture of Q that the CLG 
and the CTTP belong to the same universality class pj§ . Analogous to the above 
presented analysis we derive the mean-field critical behavior of the order parameter 
of the CTTP with random neighbor hopping. 

In the following we denote the densities of sites with p a (double occupied and 
active), pi (single occupied and inactive), and p c (empty). Normalization requires 
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Table 2. The configuration of aCTTP lattice before (s,ti,t 2 ) and after (s',^,^) 
a particle hopping. Only the source lattice site (s) and its two targets sites (ti and 
t2) where the two particles may be moved are shown. An denotes the change of 
the number of active sites due to the particle hopping and p is the corresponding 
probability of the configuration (s, ti,t2) if one neglects spatial correlations. 
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Pc + Pi + Pa = 1 and the particle conservation is reflected by the equation p\ + 2p a = p 
where the control parameter p denotes again the density of particles on a D- 
dimensional lattice, i.e., p — N/L D . The probability that a given lattice site s is 
active is therefore p a . In this case the two active particles are tried to transfer to 
two randomly chosen lattice sites ti and ti- In the case that both sites are empty 
the two particles are moved to the empty sites and the number of active sites is 
decreased by one (An = —1). The probability for this process is p a pl- All other 
possible configurations and the corresponding probabilities are listed in table|. The 
probabilities that the number of active particles are changed by An are thus given by 

Pa„=-i = Pa [Po + 2p c p a ] 

Pa„=o = Pa [2 PcPi + 2 PiPa + pi] ^ 
PAn=l = PaPf- 



The steady state condition £?[An] = leads to the equations 

Pa-0 V -l + 2p-4p a + pl = 0. (17) 

Again the first equation corresponds to the absorbing state and the second equation 
yields the order parameter as a function of the particle density 

p a = 2± ^5-2p. (18) 

Here, the + solution can be neglected (p a > 1) and the — solution describes the order 
parameter behavior above the critical density p c = 1/2. Close to this critical point 
the order parameter scales in leading order as 

Pa = \s P + 0(Sp 2 ), (19) 

i.e., the non-universal metric factor of the CTTP is ci = 1/4 and the critical exponent 
is in agreement with the CLG model (3 = 1. 

Similar to the CLG model we now apply an external field which activates single 
occupied sites. The probability that the external field h acts to a given site is p\ h 
and one tries to transfer this particle to a randomly chosen lattice site. In the case 
that the activated particle is moved to an empty lattice site the number of active site 
is unchanged by this field induced process (An = 0). The number of active sites is 
increased only if the particle is moved to a single occupied lattice site (An = +1). The 
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probability for this process is pf h. In order to incorporate the external field into the 
dynamics one has to modify p An= i accordingly and the steady state condition yields 

p a (-1 + 2p - 4p a + pi) + h (p - 2 Pa f = 0. (20) 

At the critical density p c ~ 1/2 the order parameter scales with the external field 
according to 

Pa = \ h 1/2 + 0(h), (21) 

i.e., the critical exponent is again a = 2 and the non-universal metric factor of the 
CTTP is given by c 2 = 1/16. 

In order to obtain the universal scaling function R we set p = p c + p c 8p and 
transform equation (EOI) to 



Focusing to the critical scaling behavior (h — > 0, Sp — * 0, p a — > where again p^/yh as 
well as Sp/y/h is kept constant) we can neglect all irrelevant terms and get in leading 
order 



Pa / Sp p a 



1 



Vh \Vh VhJ 4 
This equation can be easily solved and one gets 



4 ^7= + 7 = °- ( 23 ) 




Pa {8p,h) = -f ± W — + -f (24) 



where the — sign can be neglected since it yields negative values of the order parameter. 
Using the non-universal metric factor c\ — 1/4 and ci = 1/16 we get eventually again 
equation (|l4|), i.e. both models the CLG as well as the CTTP are characterized by the 
same universal function R(x, y) in the mean-field solution. Furthermore, the obtained 
universal function R agrees with that of the mean-field solution of directed percolation 
(see for instance |1^|), i.e., although the CLG and CTTP differ from the directed 
percolation scaling behavior in low dimensions they coincide on the mean-field level. 



4. Numerical simulations 



In the following we compare our results with those obtained from numerical simu- 
lations. The upper critical dimension of the universality class of absorbing phase 
transitions with a conserved field is D c = 40. Thus we compare our results with 
the scaling behavior of the CLG in D = 5 §, the CTTP in D = 5 and D = 6 Q, 
as well as with the scaling behavior of a two-dimensional CLG on a square lattice 
where active particles are moved to randomly chosen lattice sites | [L0| . In all models 
the order parameter is determined as a function of the control parameter for various 
fields and the data are rescaled according to equation (|l^). Varying the non- universal 
metric factors we observe a data-collapse with the universal function R(x,l). The 
corresponding curves are presented in figure^. As one can see, all numerically ob- 
tained curves fits well with the derived universal function. Furthermore the perfect 
data collapse of the curves for different dimensions, as well as for a mean-field model 
clearly confirms that four is the upper critical dimension. 
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Figure 2. The mean-field universal function R(x, 1) [see equation (jl5|)] of the 
universality of absorbing phase transitions with a conserved fielcLThe numerical 
data of the five and six dimensional models are obtained from fel, hl| . Additionally 
we plot the data of a (mean-field like) CLG model with random neighbor hopping 
on a square lattice (z = 4 next neighbors) which was introduced in M. At least 
four different field values are plotted for each model. 



Notice that the mean-field behavior of the CTTP order parameter was recently 
considered in [T^ |. Using a cluster approximation method the authors obtained 
equation that describes the zero-field behavior of the order parameter. 
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